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ABSTRACT 

We describe the search for Lyman-break galaxies (LBGs) near the sub-millimeter bright starburst 
galaxy HFLS3 at z=6.34 and a study on the environment of this massive galaxy during the end of 
reionization. We performed two independent selections of LBGs on images obtained with the Gran 
Telescopio Canarias (GTC) and the Hubble Space Telescope (HST) by combining non-detections in 
bands blueward of the Lyman-break and color selection. A total of 10 objects fulfilling the LBG 
selection criteria at z>5.5 were selected over the 4.54 and 55.5 arcmin 2 covered by our HST and GTC 
images, respectively. The photometric redshift, UV luminosity, and the star-formation rate of these 
sources were estimated with models of their spectral energy distribution. These z~6 candidates have 
physical properties and number densities in agreement with previous results. The UV luminosity 
function at z~6 and a Voronoi tessellation analysis of this field shows no strong evidence for an over¬ 
density of relatively bright objects (ihfio5ve< 25.9) associated with HFLS3. However, the over-density 
parameter deduced from this field and the surface density of objects can not excluded definitively the 
LBG over-density hypothesis. Moreover we identified three faint objects at less than three arcseconds 
from HFLS3 with color consistent with those expected for z~6 galaxies. Deeper data are needed to 
confirm their redshifts and to study their association with HFLS3 and the galaxy merger that may 
be responsible for the massive starburst. 

Subject headings: galaxies: formation - galaxies: high-redshift - galaxies: luminosity function, mass 
function - galaxies: starburst 
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INTRODUCTION 

One of the most important questions of modern as¬ 
tronomy is undoubtedly the formation and evolution of 
the first luminous objects in our Universe. During the 
last decade, the arrival of new facilities with capabili¬ 
ties to push even further the boundaries of our Universe, 
such as for example Herschel Space Observatory (|Pil- 


brattetal. 2012), the HST/WFC3 (|Robberto eta 


2002]), and vLT/HAWK-I (Pirard et al. 2004), con¬ 


tributed to considerable advances m our understanding 
of the first billion years. The number of z >6 galax¬ 
ies presently confirmed by spectroscopy h as stro ngly in¬ 
creased over the last few years (e.g. Schenker et al. 
2012, Vanzella et al. |201l|) allowing to better constrain 
their physical properties in terms of the star -format ion 
rate ( S FR), stellar mass, and reddening (e.g. de Barros 
et al. ||2012). The evolution of such galaxies is now rel¬ 
atively well-constrained and supported by spectros copic 
observations out to z~6 (e.g. Le Fevre et al. , and 

photometric studies are now slowly starting t o give ro¬ 
bust co nstraints of galaxies out to z ~10 (e.g. 
et al. 20T§. 


Bouwens 


The detection of large numbers of dusty, massive star 
burst galaxies at z ^2 was a surprise when they were 
f irst id entified in 1997 (Smail et al. ||1997l Hughes et al. 


1998 for a review see Casey et al. 2014). Theoreti¬ 


callythe existence of dusty, massive star bursts at such 
early epochs is difficult to explain and thorough observa¬ 
tional constraints on their properties provi de a str ingent 
test of the galaxy formation models (e.g. Baugh et al 


2005). The HerMES survey (Oliver et al. 2012| ) ha's 


id entified a population of dusty star burst galaxies at z > 


4 (Dowell et al. ||2014) that are n ot predicted by today’ s 
galaxy formation paradigms (e.g. Hayward et al. |2012) 


The efficient identification and detailed study of these 
rare z > 4 star-bursting galaxies is important for under¬ 
standing t heir progenit or populations, drivers of cold gas 
accretion (Carelli et al. 2010), and their descendants. It 
has been shown by previous studies that several massive 
starburst at high -z could be proto-cluster members, but 
none of them has been identified at z >6 (Daddi et al. 
2009} | Capak et al. ||MTT ) 


Recently, Riechers et a. 


(2013) have discovered an ex¬ 


treme sub-millimeter starburst at z =6.34, called HFLS3 


hereafter, ov er 21 deg 2 of the Herschel /SPIRE ([Griffin 
et al. 2010) data from the HerMES survey. This object 
is lensed by a two-component galaxy s ystem at z ^2.1 
involving a magnification of 2.2 ±0.3 (Cooray et al. 
2014). Its physical properties based on the rest-frame 
UV emission results in a SFR of 1320 M 0 yr _1 and dust 
and stellar masses of 3 xlO 8 M 0 and 5 xlO 10 M 0 , re¬ 
spectively, making this source one of the most massive 
starburst presently known during the epoch of reioniza¬ 
tion (EoR). 

A recent paper aiming to detect sub-millimeter emis¬ 
sion in the vicini ty of HFLS3 with SCUBA-2 data (Rob¬ 


son et al. 2014) found no evidence in favor of an over¬ 
density of dusty galaxies associated with HFLS3. In 
this paper, we present results of an optical and near- 
IR analysis of the environment of HFLS3 aiming to re¬ 
cover Lyman-break galaxies at z~6. We combine two 
datasets: ground-based images covering a wide field of 
view (to select the brightest objects at z~6) and deeper 


Hubble/ WFC3 and ACS data with a smaller field of view 
(to detect fainter objects). 

In Section [l] we present our new data used for the 
analysis and their reduction procedure. In Section [2j 
we explain in detail the method that was used to se¬ 
lect z~6 objects and we present the candidates coming 
from the two ground and space-based datasets. Their 
photometric properties are presented in Section [3] In 
Section [4] we discuss the presence of an over-density of 
LBGs in the HFLS3 field. The concordance cosmology 
is adopted throughout this paper, with Ua= 0.7, U m =0.3 
and Hq =70 km s" 1 Mpc" 1 . All magn itudes are quoted 
in the AB system ( Oke and Gunn [1983 ). 


1. OBSERVATIONS AND DATA REDUCTION 

In this Section, we present the ground-based and space- 
based data we used to perform the search for z~6 objects 
around HFLS3. GTC optical and Spitzer dat a have al¬ 


ready been presented in Riechers et al. (2013). We have 


re-reduced the GTC data to improve their image qualities 
and we include a newly acquired K s -band image taken 
with LIRIS on the the 4.2m William Herschel Telescope 
(WHT). In the following we summarize image properties 
and the data reduction procedure. A summary of the 
data used in the current study is listed on Table [l] 

1.1. Ground-based data 

In order to select bright objects at z~6 and to esta b 


lish the bright-end of the UV luminosity function (Lacey 


et al. 2011), we performed a first search for LBGs in 


ground-based images taken with O SIRIS on GTC ([Al 
varez et al. 1998, |Cepa et ah 


. 2003]) and LIRIS installed 

on WHT (| Wall et al.' ||1987[ |Manchado et al. |[l998|) . 

OSIRIS data were acquired between June 29 and Au¬ 
gust 3 2011, as part of the GTC2-10ITP program (P.I.: 
I. Perez-Fournon), and used g ', r', i'and z' broad-band 
filters. The field of view covered by OSIRIS is 7.8’ x 
7.8’ and the total exposure times in each band are 2.7 
ksec in g' and r', 21.6 ksec in i' and 18,7 ksec in z '. Each 
frame was reduced individually following standard reduc¬ 
tion procedure in IRAfj 2 ^] (bias subtraction, flat fielding, 
sky subtr action) and registration and co mbination usin g 
SCAMP QBertin E. | |2006) and SWARP ( |Bertin |[2010| ). 
To reduce the sky background at long wavelengths we 
applied for the i' and z' band images, a double sky sub¬ 
traction using thn IRAF pac kag e XDIMSUM, with steps 
explained in Richard et al. (2006). The astrometric and 


and LIRIS insta 


photometric calibrations wer e performed using the 9 th 


release of the SDSS catalog (Ahn et al. 2012). More¬ 


over to produce the best z— band quality image, we ap 
plied the ’’best seeing stacks” method described in |Gwyn| 
(2012), selecting only the best frames in terms of seeing. 
We measu red the FWHM on each frame using PSFEx 
(|Bertin 20lT|, and we kept 86.4% of our original dataset 
beyond which the FWHM increase without a significant 
evolution of the image depth. The seeing of the final im¬ 
ages is 0.89”, 0.83”, 0.83” and 0.89” respectively in g\ 
r', i' and z'. The depth was computed using empty 1.4” 
radius apertures all over the field. 

We added to our previous dataset a new near-IR, K s 
image acquired on June 5 2014 with the LIRIS instru- 

26 http://iraf.noao.edu/ 
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ment (4.27’x4.27’ field of view, 0.25” size pixels) in¬ 
stalled at the Cassegrain focus of the WHT, as part of 
a new WHT Large Program on HerMES high-redshift 
galaxies (P.I.: I. Perez-Fournon). The total integra¬ 
tion time was 3.6 ksec, the observing conditions were 
photometric and the seeing was 0.5”. The reduction of 
the image was performed using the I AC’s IRAF lirisdr 
taskpjand the astrometry calibration was carried out us¬ 
ing the Graphical Astronomy and Image Analysis Tool 
(GAIA ) inclu ded in the Starlink astronomical software 
packageThe K s image was matched to SDSS DR10 
stars and is estimated to be accurate to 0.1” (rms of the 
fit). The photometric calibration of the image was done 
with 2MASS stars. The 5-sigma depth of this image in 
1.4” radius apertures is 23.14. 


1.2. Space-based data 

The data reduction of the ACS and WFC3 data we 
used in this study are described in detail in Coor ay et| 
al. (2014). We summarize hereafter the principal steps. 
These data are part of an HST Cycle 21 program (GO 
13405, PI: A. Cooray) and used 6 orbits to produce 
images of the HFLS3 environment in F625W, F814W, 
F105W, F125W and F160W (Table 0. We used the 
IRAF.STSDAS pipeline and the CALWFC3 tool for flat- 
fieldi ng and cosmic-ravs rej ections and ASTRODRIZ- 


ZLE (Fruchter et al. 2010) to combine individual ex¬ 


posure and create 0.06^/pixel WFC3 images. The ACS 
data were flat-fielded, charge transfer efficiency corrected 
and cosmic rays rejected using the CAL ACS pipeline 
(version 2012.2), and then combined with ASTRODRIZ- 
ZLE to produce 0.03”/pixel images. Astrometric calibra¬ 
tion was performed using SDSS ( 9 th release). The WFC3 
and ACS data were matched independently following the 
classical IRAF procedure. The depth was computed us¬ 
ing empty apertures all over the field (0.25” radius for 
ACS and 0.50” radius for WFC3 images). 

In order to extend the wavelength coverage of our sur 


vey, w e also used data obtained with IRAC (Fazio et al 
2004) on board the Spitzer Space Telescope as part of a 
DDT program (ID: 80240 - PI: J. Vieira) on 2012 March 
21. The on-target observations consisted of 38 frames 
with an integration time of 100 sec in each of the 2 chan¬ 
nels centered at 3.6 /am (chi) and 4.5 /am (ch2). We used 
corrected Basic Calibrated Data (cBCD) frames which 
are already corrected by pipeline for various artifacts, 
such as multiplexer bleed and pulldown. These images, 
together with associated mask and uncertainty images, 
were processed, drizzled (with a PIXFRAC of 0.65), and 
combined using the standard pipeline MOPEX. The final 
mosaics images have a pixel size of 0.6”, roughly a half of 
the IRAC native pixel scale. The mosaic has a 3cr point 
source sensitivity of 0.384 /iJy and 0.412 /iJy in 3.6/im 
and 4.5/im respectively. 


1.3. Source catalogs 


We used SExtractor 2.18.4 (Bertin & Arnouts 1996) 
to build catalogs in double image mode using a y 7 im¬ 
age including all individual picture, as detection pic- 



Fig. 1.— Wavelength coverage and filters transmission. 


TABLE 1 

Photometric properties of the imaging data 


Filter 

V// 

AA 

^exp 

m(5<r) 

Instrument 


[nm] 

[nm] 

[ks] 



(i) 

(2) 

(3) 

(4) 

(5) 

(6) 

F 

481.5 

153 

2.7 

26.91 a 

OSIRIS/GTC 

r’ 

641.0 

176 

2.7 

26.80 a 

OSIRIS/GTC 

V 

770.5 

151 

21.6 

26.60 a 

OSIRIS/GTC 

z ’ 

969.5 

261 

16.3 

25.87 a 

OSIRIS/GTC 

F625VF 

629.6 

98 

1.6 

26.08 6 

ACS/HST 

F814VF 

811.5 

166 

2.3 

26.99 b 

ACS/HST 

F105VF 

1055.1 

265 

9.9 

25.90 c 

WFC3/HST 

F125VF 

1248.6 

443 

4.4 

26.30° 

WFC3/HST 

F160VF 

1536.9 

268 

2.8 

26.00° 

WFC3/HST 

Ks 

2150.0 

320 

3.6 

23.14 a 

LIRIS/WHT 

3.6//m 

3575.0 

776 

3.8 

24.38 d 

IRAC/Spitzer 

4.5/xm 

4528.0 

1060 

3.8 

24.30 d 

IRAC/Spitzer 


Note. — (1) filter identification, (2) filter central wavelength, 
(3) filter width, (4) exposure time, (5) 5cr AB magnitude and (6) 
instrument and telescope. 

a point-source depth measured in a 1.4” radius aperture 
b point-source depth measured in a 0.25” radius aperture 
c point-source depth measured in a 0.50” radius aperture 
d point-source aperture corrected depth measured in a 1.9” radius 
aperture 


ture (Szalay et al. 1999). The extraction param¬ 
eters we used have been chosen in order to extract 
faint objects. Therefore the DETECT_MINAREA pa¬ 
rameter was fixed to 4.0 pixels according to the see¬ 
ing of our images, and to limit spurious detection s we 
choose DETEC T_THRESH =1.7cr (e.g. see 
et al. p 
M 


Muzzin 

20131 DEBLEND_NTHRESH=16 and DE- 
ELEN D_MlNUDNT=0. 00002 flScoville et al. |2007|). We 
measured the photometry in apertures with radius de¬ 
fined by 2xFWHM on PSF-matched data on the GTC 
data and 0.4” radius aperture on HST images, as well 
as with SExtractor MAG_AUTO. We found that these 
two types of magnitude are in good agreement for all 
point-source like objects in our catalogs. 


27 www.iac.es/galeria/jap/lirisdr/LIRIS_DATA_REDUCTION.html 

28 http://star-www.dur.ac.uk/ pdraper/gaia/gaia.html 

29 http://starlink.jach.hawaii.edu/starlink 


2. SELECTION OF Z~6 LYMAN BREAK GALAXIES 


The LBG technique (Steidel et al. 11996) hasjoeen ex¬ 


tensively used to select galaxies at high-z (e.g. Castellano 
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et al. [2010 Trenti et al. J|2012 , Laporte et al. 1 2014). It 


combines two different criteria: non-detection / detection 
and the color sel ection. We applied this m ethod on SEx- 
tractor catalogs ( Bertin fe Arnouts 11996 ) built using de¬ 
tection parameters as a function of the dataset we used 
(see below for details). Regarding the depth and spa¬ 
tial resolution of each dataset, we performed independent 
searches in each survey: from our ground-based images 
we focussed on the brightest objects and from the HST 
data we investigated faintest sources at the redshift of 
HFLS3. 


2.1. Selection from GTC data 

The GTC/OSIRIS imaging data we used to search 
for LBG at the redshift of HFLS3 covered a wave¬ 
length range between 328.5 and 1230.5 nm. The non¬ 
detection/detection criteria we applied to our catalogs 
to select objects were the following: 


m g > > m(2a) g ' 

= 27.9 

(1) 

m F > m(2cr) r f 

= 27.8 

(2) 

m z / < m(5cr)^/ 

= 25.9 

(3) 


The lack of deep near-IR data over the OSIRIS field-of- 
view does not allow us to produce a color-color diagram 
with a clear selection window. 

We defined a color window using the i' — z' color evo¬ 


lutio n of a set of templat es (Bruzual & Chariot 2003 
and |Polletta et al. ||2007|) as a function of the redshift 
(figurd2l rrq-m^ >0.0). These criteria limit the contam¬ 
ination by most of the low -z interlopers except for the 
elliptical galaxies at z <3. In order to avoid the contam¬ 
ination by such interlopers, we estimated the size of the 
break between optical and NIR beyond which the selec¬ 
tion of mid ,-z interlopers is unlikely (figure [3]). The break 
required by the evolution of templates mentioned above 
(m r -m z >2.5 mag) is consistent with the non-detection 
criteria that we applied, and thus by combining the pre¬ 
vious i' — z' color criterion (nq- m z >0.9 mag) with the 
non-detection/detection criteria listed above, we elimi¬ 
nate a large number of mid -2 interlopers. Therefore, we 
can summarize the color selection we defined to select 
z~6 objects on the GTC data as follows: 


nq - m z > 0.9 
m r - m z > 2.5 

The selection window is shown on Figure [4j 
After visual inspection only two objects satisfy the se¬ 
lection criteria defined above with m z <25.1. These two 
candidates are displayed as zl_GTC and z2_GTC in Fig¬ 
ure [U and Tabled 


2.2. Selection from HST data 

The second sample has been built using HST data cov¬ 
ering wavelength from 580.6 to 1670.9 nm. In order to 
select z~6 objects and regarding the set of filters we have, 
we defined the following detection/non-detection criteria: 

mF625W > m(2cr)F625W U m F i 05W < m(ba) F1 o 5W (4) 
or in terms of color : 


mF625W — r MFlG5W > 1-2 . (5) 

The wavelength coverage of the HST data is better in the 
near-IR domain compared to the GTC survey. There¬ 
fore we can define a color-window to select z~6 objects 



o - XJ -^ ^ I - * -1-i- 

0 2 4 6 8 


Z 

Fig. 2. — i > -z > color cut defined for the z~6 LBG selection using 
GTC filters and color evo lutio n of several template s of galaxies 
(|Bruzual & Chariot ||2003| and |Polletta et al. ||2007|). The color- 
criteria is defined by the black line, and shows tfiat in that case 
the majority of the contaminants are elliptical galaxies at mid-zs. 



Fig. 3. — r'-z* color window defined for the z~6 LBG selection 
using GTC filters. 


following the standard method described in section |2.1| 
Figure [6] displays this selection window defined by: 

r rriF625W — mF8i4W > 2.0 ( 6 ) 

r m F 8,14W - rnFlOSW < 4.0 (7) 

TTlF625W — JTIF814W > 0.57 X (m F si4W ~ mF105W ) + 0.66 

( 8 ) 

The use of broad-band filters as well as color selection 
imply that several of the z~6 candidates will be lost dur¬ 
ing the selection process. In order to estimate the frac¬ 
tion of those objects, we computed the incompleteness 
of the selection function described in this paper. Us¬ 
ing templates from Coleman et al. (1980), Kinney et al. 
(1996), Poletta el al. (2007) and Silva et al. (1998) and 
transmission of filters used in this survey, we simulated 
100 000 objects in the redshift range 4.0 < 2 < 8.0 per 
bin of 0.25 magnitude with Hip 10 svf ranging from 19.0 
to 29.0. The selection by the color-criteria defined above 
shows that ^74% of objects at z~6 .3 are selected. 

As previously we used SExtr actor 2.18.4 in doub le im¬ 
age mode using a y 2 image (Szalay et al. 1999) made 
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TABLE 2 

Photometry of the two candidates selected from the GTC/OSIRIS data 


ID 

RA 

DEC 

m g 

m r 

m i 

m 2 

m K.s 

m 3.6 

m 4 .5 

zl_GTC 

17:06:40.6 

+58:47:49.6 

>27.9 

>27.8 

27.92 

24.96 

23.08 

23.24 

24.73 






±0.77 

±0.10 

±0.22 

±0.08 

±0.22 

z2_GTC 

17:06:41.310 

+58:47:17.28 

>27.9 

>27.8 

27.48 

25.05 

23.62 

22.91 

23.00 






±0.51 

±0.11 

±0.36 

±0.20 

±0.20 

Note. — 

All GTC-magnitudes are MAG_AUTO magnitude from SExtractor. Upper limits are 2cr depth. 

Error bars ; 

are computed within empty 1.4” 

radius apertures around the object 

on the GTC data. IRAC 

magnitudes 

are measured in a 1.9” radius aperture and error 

bars are 

computed using empty 1.9’ 

’ radius 


aperture around the object. 



J_ L _i_I_i_ 

0 2 4 


i-z 


Fig. 4. — Color criteria used to selec t z >5.5 objects using tem - 
plates from |Bruzual & Chariot | (|2003|) and |Polletta et al. | (|2007| . 
Position of the two OTU candidates are displayed by magenta dots. 
We over-plotted the expected color of M, L and T dwarfs in grey 
computed from 225 stellar spectra (see references in Sec. |2.4| ). 
Black dots show the colors of objects with redshift ranging from 
5.5 to 6.5 per bin of 0.1 


with the WFC3 images as detection picture to produce 
WFC3 catalogs, because LBGs should be detected 
at these wavelengths. As in the ” ground-based” selec¬ 
tion, we used SExtractor parameters defined to maxi- 
mize t he select ion of faint objects. We then used TOP- 
CAT ([Taylor |[2005|) to match ACS single image mode 
catalogs with the VVFC3 catalogs allowing a search ra¬ 
dius fixed to 2xFWHM of the worst seeing image. The 
total magnitude were obtained from aperture correction 
using F160W MAG _ AUTO and follow ing the method de¬ 
scribed in Finkelstein et al. (2013). Each catalog in¬ 
cluded ~1600 detections. Colors were measured on psf- 
matched data, and 14 objects follow the criteria defined 
above. Among these sources, one of our GTC candidates 
is covered by the HST field of view and is included in 
this new sample (ID: z2_GTC), showing that our two 
selection functions are well defined to select this kind 
of object. The other GTC candidate (zl_GTC) is not 
detected on the F625W images and well detected on the 
F814W (rnF8i4W= 26.80T0.18) confirming the break be¬ 
tween the %' and 2 /-band of our GTC survey. However, 
HFLS3 is not included in this sample, because SExtrac¬ 
tor failed to extract it properly, it is blended with the 2 


nearby z ~ 2.1 galax ies acting as a foreground lens (see 


Cooray et al. 


2014 for more details). Five objects have 


been removed from this sample because they are detected 
on t he d eep g' band data from the GTC survey (see sec¬ 
tion |2.4| for details). Therefore the final sample is com¬ 
posed of nine objects and is presented on Table [3] and 
Figure [7[ 

2.3. Detection in LIRIS and IRAC data 

We increased the number of SED constraints in the 
near-IR domain by using data from LIRIS and IRAC. 
Due to the depth of our K s image, only the two GTC- 
objects are clearly detected in the near-IR (including the 
candidate in common in the two samples). These re¬ 
sults are not surprising regarding the brightness of these 
sources, but it could have been a way to remove mid -2 in¬ 
terlopers from the HST sample. Indeed previous papers 
have shown that mid .-z interlopers display a very high-z 
LBG like SED up to 1.6/im, but they are very red at 
larger wave lengths corresponding to IRAC (e.g. Boone 
et al. 2011) The depth of IRAC data (m 3 . 6? 4.5(5cr)^24.3) 
is not completely sufficient to add robust constraints on 
the SEDs of all our candidates. Among the two different 
samples, two objects are clearly detected in the IRAC 
data. We extracted the photometry of z2_GTC using 
^ For both objects, we measured the photometry 


in a 1.9” radi us aperture, and used the correction factor 
computed by Mauduit et al. (2012). The error bar was 
computed using empty 1.9” radius apertures around the 
object. 


2.4. Contaminants 

The most common sources of contamination of high- 
z samples are SNe, AGN, low-mass stars, photometric 
scatter, transient objects, spurious sources, extremely 
red galaxies, among others. All our candidates are de¬ 
tected in at least two bands limiting the contaminations 
by spurious sources, moreover these observations spread 
over several months limiting the detection of transient 
sources. However, we noted that our ample is not free of 
AGN and SN contamination given the timescale during 
which our observations have been carried out. Regard¬ 
ing the field of view covered by these two surveys and 
the statistical number of supernovae expected in ^ldeg 2 , 
the probability of being contaminated by a supernova is 
relatively small. Low-mass stars, such as M, L and T 
dwarfs, could display the same colors properties as those 

30 www.mpe.mpg.de/ erwin/code/imfit/index.html 
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Fig. 5.— Thumbnail images of the two candidates selected in the 7.8’x7.8’ OSIRIS field of view. The size of each stamps is 12” xl2”. 
The position of the candidate is displayed by a white circle of 1.4” radius aperture. 

into account. 

Therefore the majority of contaminants that can enter 
into our sample involves the extremely red mid r z galax¬ 
ies. We estimated the contamination rate of our sample 
by mid rz interlopers by using the colors distribution of 
mid -z objects over the fields of view covered by our two 

2010). We followed a 3 steps 
(1) we selected all ob- 



datasets (e.g. Oesch et al. 


method summarized herea: 


ter 


Fig. 6.— F625W-F814W-F105W color window defined for the 
z^6 LBG selection using HST filter s and color evo l ution of se veral 
templates of galaxie s given by|Bruzual & Chariot |(|2003|) and|Pol-| 
|letta et al. | (2007|). The color-criteria tor z >5.5 are defined by 
the black box. We over-plotted the expected color of M, L and T 
dwarfs in gr ey computed from 225 stellar spectra (see references 
in Sec. \2A\ . Black dots show the colors of objects with redshift 
ranging from 5.5 to 6.5 per bin of 0.1 

expected for high-z galaxies. We computed the expected 
colors of these galaxies through all the fil ters we used in 
this st u dy from 225 stellar spectra 


Burgasser et 


al. 


2006b , Burg asser et al. ||2004[ |Burgasser et al. | 2008 


al. 2004 

Kirkpatrick et al. 

2010 ; 


Reid et a 


2006 

Siegler et a 

1. 2007 

Chiu et al. 

2006, 

^ooper et al. 

2007 


et al. ||2007[ |Burgasser et al. ||2006c|). Fig. [6] and [4 

display the expected colors of t hese stars that are consis¬ 


tent with previous s tudies (e.g. Finkelstein et al.~||2014 


Willott et al. |20l3 ) and show that stellar contamination 


is non-negligible is both samples. Therefore to limit the 
selection of low-mass stars, expected to be unresolved on 
our dataset, we used SExtractor FLUX_RADIUS enclos¬ 
ing 50% of the total flux and the SExtractor stellarity 
parameter. Both indicators demonstrate that two ob¬ 
jects among our samples are likely unresolved, namely 
zl_HST and z3_HST with a stellarity parameter of > 0.7 
and a size computed from the half light-radius compara¬ 
ble with the FWHM of the image, taking uncertainties 


jects that are detected at more than 2a in all bands, (2) 
then we matched the luminosity range to the luminosity 
covered by our two samples and added the corresponding 
uncertainties on the photometry, (3) finally we applied 
the selection criteria we used to build our high -z can¬ 
didates samples and all the selected objects are mid-2: 
interlopers. 3 objects are identified over the 4.5 arcmin 2 
covered by HST data and 2 using our GTC data leading 
to ^ 7-36% and >20% contamination rate including cos¬ 
mic variance respec tively. It has been shown by several 


recen t studies (e.g. Hayes et al. 2012, Laporte et al. 


2015) that the best way to remove such objects is to cover 


the short wavelength with deeper imaging data, regard¬ 
ing that the contamination rate by such extreme mid-2: 
interlopers is rather uncertain. We used our deep GTC 
g'-band data to confirm the real non-detection of HST 
candidates. We find that one of our targets is clearly de¬ 
tected (m g ^4 a) and can not, therefore, be at z~6. Four 
others candidates are faint on that image (m p / ~2a) and 
thus cannot also be at such high- z, but further data are 
necessary to confirm their optical emission. 

After removing these mid-2: interlopers, the final HST 
sample is composed of nine objects with mpios ranging 
from 24.5 to 26.5 (figurej8l). The postage stamps are dis¬ 
played on Figure [7] and Table [3] presents the photometry 
of this sample. 

We note here that our ACS images are not deep enough 
to exclude definitively the low-2: hypothesis for all our 
2:^6 galaxy candidates. Indeed, the 2 a depth of the 
F625W image used here (27.1 AB) combined with the 
F814W depth (27.99 AB) does not allow us to apply 
completely the color-color criteria we imposed (F625W 
- F814W >2.0). For most of them, the F625W-F814W 
break is not high enough and even if we used the la de¬ 
tection as upper limit in F625W we never reach a break 
of 2 magnitu des i n thes e two filters (as in previous stud¬ 
ies, e.g. Monna et al. 2014). This is why we primarily 


use here the deeper ground-based GTC g-band data to 
identify the interlopers. 
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F625W 

ACS/HST 


F814W 

ACS/HST 


F105W 

WFC3/HST 


F125W 

WFC3/HST 


F 160 W 

WFC3/HST 


Ks 

LIRIS/WHT 


3.6|jm 

IRAC/Spitzer 


4.5|jm 

IRAC/Spitzer 



zl HST 


z2 GTC 


z2 HST 


z3 HST 
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z5 HST 


z6 HST 


z7 HST 


z8 HST 


z9 HST 


Fig. 7.— Thumbnail stamps image of the nine candidates selected in the WFC3 field of view. The size the stamps are 7.5” x7.5” in all 
the bands. The position of the candidate is displayed by a white-circle of 0.8” radius aperture . 


TABLE 3 

Photometry of the candidates selected in the HST data 


ID 


RA 


DEC 


m F 625 


m F8i4 m^i05 125 


m 3.6 


m 4 .5 


zl_HST 17:06:41.310 +58:47:17.28 >27.08 


z2_HST 

17:06:40.737 

+58:46:58.52 

>27.08 

z3_HST 

17:06:57.157 

+58:46:31.85 

>27.08 

z4_HST 

17:06:40.411 

+58:47:01.39 

>27.08 

z5_HST 

17:06:55.049 

+58:46:39.05 

>27.08 

z6_HST 

17:06:40.212 

+58:45:57.16 

>27.08 

z7_HST 

17:06:45.315 

+58:46:39.78 

>27.08 

z8_HST 

17:06:45.289 

+58:45:20.86 

>27.08 

z9_HST 

17:06:49.819 

+58:46:58.26 

>27.08 


26.42 

24.52 

24.09 

23.90 

22.91 

23.00 

± 0.13 

± 0.07 

± 0.03 

± 0.03 

±0.20 

±0.20 

>27.99 

25.08 

25.37 

25.08 

>25.38 

>25.30 


± 0.11 

± 0.10 

± 0.09 



27.01 

25.34 

25.07 

25.09 

>25.38 

>25.30 

± 0.23 

± 0.13 

± 0.07 

± 0.09 



27.91 

25.81 

25.53 

25.53 

>25.38 

>25.30 

± 0.47 

± 0.20 

± 0.11 

± 0.14 



>27.99 

25.87 

25.61 

25.60 

>25.38 

>25.30 


± 0.21 

± 0.12 

± 0.15 



>27.99 

25.92 

25.64 

25.81 

>25.38 

>25.30 


± 0.23 

± 0.12 

± 0.18 



26.73 

25.94 

25.67 

25.68 

>25.38 

>25.30 

± 0.15 

± 0.23 

± 0.13 

± 0.16 



27.67 

26.11 

25.80 

26.01 

>25.38 

>25.30 

± 0.37 

± 0.27 

± 0.14 

± 0.22 



26.86 

26.47 

26.45 

26.43 

>25.38 

>25.30 

± 0.18 

± 0.37 

± 0.25 

±0.32 




Note. — All HST-magnitudes are MAG.AUTO magnitude from SExtractor. Upper limits are 2a depth. Error 
bars are computed within empty 0,125” radius apertures around the object on the ACS data and 0.25” radius on 
WFC3 data. IRAC magnitudes are measured in a 1.9” radius aperture and error bars are computed using empty 
1.9” radius aperture around the object. 
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2.5. Expected number of objects 

The effective surface covered by our OSIRIS survey is 
55.5 arcmin 2 and has been computed by masking all the 
bright objects on our y 2 detection picture. We computed 
the expected number of sources in that field of view by 
using the shape of the l uminosity function published by 


Bouwens et al. 
lished at 


Tnc 


(2012), which is relatively well estab- 
uding cosmic variance and uncertain¬ 
ties on the Schechter parameters, the expected number 
of objects in the range of redshift 5.5< z <6.5 and with 
m z <25.9, is 1.3+o.g, showing that our sample built from 
GTC data is in excellent agreement with the expecta¬ 
tion. For the HST-survey covering an effective surface of 
4.54 arcmin 2 , the expected number of objects detected 
at most than 5cr with WFC3 bands in the same area is 
6.10^4.1 (including cosmic variance). Therefore our se¬ 
lection is in relatively good agreement with the expecta¬ 
tions, without any evidence for a significant overdensity. 

3. PHOTOMETRIC PROPERTIES OF THE Z~6 
CANDIDATES 

This Section presents the principal photometric prop¬ 
erties of the selected objects using a SED-fitting ap¬ 
proach. We discuss in the following the photometric red- 
shift, UV luminosity, SFR, and the reddening for our two 
samples. 

We use d the Version 12.2 of Hyperz^] (Bolzonella et 
al. [20 00) with the standard templates library including 
Bruzual & Chariot (2003), Coleman et al. (1980), Kin¬ 
ney et al. (1996), Poletta el al. (2007) and Silva et al. 
(1998) templates. The redshift space parameters ranges 
from z ^0 to 8 and A v from 0.0 to 3.0 mag. The best-fit 
is always found at 2 >5.3 with a reasonable reddening 
solution for most of the candidates (A v <1.00mag). The 
redshift probability distribution, hereafter P(z), is well 
defined around the best photo-z and no low -2 solution 
appears clearly for all our targets. For the candidate in 
common in the two samples, we estimated its properties 
in three cases: (1) combining all of the data points we 
have, (2) using only the GTC, WHT and Spitzer pho¬ 
tometry, and (3) using the HST, WHT and Spitzer con¬ 
straints. The results are similar and give a moderate 
reddening solution (A v ^1.80-2.0) and a best photo-z 
included between 5.8 and 6.1. 

As a cross-check, we also tried to fit the SEDs of all 
these dropouts assuming a low-z solution, i.e. with a 
redshift space parameters ranging from 0.0 to 3.0 and al¬ 
lowing the same reddening interval. The best-fit is found, 
in each case, with a higher y 2 and a P(z) not well de¬ 
fined in the redshift interval explored (cf. Table [4J. We 
also used a prior in luminosity to check the consistency 
of our sample with previous studies at z~6. The prior 
was defined using the parametrization of the Schechter 
function published in Bouwens et al. (2014). The pa¬ 


rameters space was defined as previously and the same 
templates library was used. Two candidates showed dif¬ 
ferent best-fit photo-Vs, but the best-fit SED template 
was always found at z > 5.5. 

The UV luminosities were computed using the best-fit 
SEDs found by Hyperz and the SFRs were computed 
from the UV luminosities u s ing th e relationship pub¬ 
lished in Kennicutt et al. (1998). We corrected for 


the dust extinction of L 1500 following the Calzetti et al. 
| (2000) method. We noticed that the I /1500 is ranging 
from ~2.5 to 29 xl0 41 ergs s -1 A -1 and the SFRs from 
^26 to 300 Mq yr -1 . These values are in good agree- 


& de Barros 

2(J10l and with properties of previous sam- 

pies at z~6-' 

7 (Zheng et al. 2009 

Curtis-Lake et al. 


The estimation of the photometric redshift generally 
implies the use of templates of nearby galaxies extrap¬ 
olated to the very high- z. We made use of the large 
number of z~6 galaxies spectroscopically confirmed to 
check the reliability of our template s library and the 
SED-f i tting method using dat a from [T oshika wa et al. 
(|2012|),|Jiang et al |(|2013) and|Willott et al. | (|2013D . We 
matched photometric catalogs published in previous ref¬ 
erences to the quality of our dataset, more especially for 
the filters used here in the wavelength range covered by 
our data (from 0.4 to 1.6 /im). We re-computed photo¬ 
metric error bars of spectroscopically confirmed galaxies 
using the average noise measured in each of our images. 
In the case of a non-detection, we used the upper limits 
computed in our data. The photometric redshift of the 
confirmed z~6 galaxies were computed with Hyperz us¬ 
ing the same templates as were used previously with the 
actual sample. Among the 50 galaxies used, 45 have a 
1<7 confidence interval for the photometric redshift that 
includes the spectroscopic redshift. The mean absolute 
dilation to the spectroscopic redshift is ^9% with a stan¬ 
dard deviation of ~16% showing that the method we 
used to estimate the photometric redshift of our candi¬ 
dates is efficient at z~6. We also used this spectroscopic 
sample to check the reliability of our color criteria, and 
demonstrated that all these spectroscopically confirmed 
galaxies would be selected using the criteria defined in 
Sec. [21 

4. ENVIRONMENT OF HFLS3 

In this Section, we describe the environment of the 
starburst galaxy HFLS3 , and discuss the possibility for 
an overdensity associated with it. 

4.1. Faint z~6 candidates over the wide field 

In order to select faintest objects at the redshift of 
HFLS3, we relaxed the detection constraints we imposed 
for our HST selection (eq. E| by requiring a detection 
>3cr in F105W and F125W (instead of the 5 <j level re¬ 
quired previously on F105W) and by using the F814W- 
F105W color criterion computed previously. Indeed by 
relaxing the detection level on the F105W image, we are 
not able to use the color criterion combining F625W and 
F814W, because of the shallower depth of the F625W. 
After visual inspection removing false detection as ex¬ 
plained above, 25 objects satisfy this new selection func¬ 
tion. In order to study if this sample of faint objects 
shows any evidence for an overdensity of galaxies asso¬ 
ciated with HFLS3, we applied the Voronoi tessellation 
method as described in|Ramella et al. Id2001) and based 

rdoi _ 


31 http://userpages.irap.omp.eu/ rpello/newhyperz/ 


on the ” triangle” C code published by bnewchuk (1996). 
The density threshold we used to distinguish background 
regions and fluctuations which are significant overdensi¬ 
ties is estimated at 1.74 object per Voronoi cell. As seen 
in figure [lOj no over-dense region is highlighted by the 
new faint objects sample. 
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TABLE 4 

Photometric properties of the z~6 galaxy candidates selected in this study 


ID 



high- z 



low- z 

l 

Physical properties 

Zphot 


A<y 

[mag] 

lcr interval 

z phot 

x 2 

A-y 

[mag] 

L 1500 

[xl0 41 erg/s/A] 

SFR 
[M© / yr] 

zl_GTC 

6.5 

0.6 

0.8 

6.3 - 6.5 

1.1 

2.2 

0.0 

28.8 

302.8 

z2_GTC 

6.1 

0.5 

1.8 

5.8 - 6.3 

1.5 

0.8 

0.8 

22.5 

236.0 

z2_HST 

6.5 

1.5 

0.2 

6.0 - 7.2 

1.4 

4.6 

0.0 

4.5 

47.0 

z3_HST 

6.2 

0.7 

0.8 

5.7 - 6.6 

1.5 

2.0 

0.0 

16.1 

169.3 

z4_HST 

6.4 

0.3 

0.8 

5.7 - 7.8 

1.5 

1.3 

0.0 

11.6 

121.8 

z5_HST 

7.1 

0.07 

0.4 

5.7 - 8.0 

1.5 

1.3 

0.0 

5.3 

56.2 

z6_HST 

7.1 

0.02 

0.2 

5.8 - 8.0 

1.5 

1.3 

0.0 

3.2 

33.2 

z7_HST 

5.6 

0.2 

0.8 

4.1 - 6.2 

0.1 

0.3 

0.4 

8.2 

86.0 

z8_HST 

6.2 

0.3 

0.8 

5.3 - 7.4 

1.5 

0.6 

0.0 

8.5 

89.5 

z9_HST 

5.2 

0.03 

0.6 

0.0 - 6.1 

0.9 

0.1 

2.2 

2.5 

26.2 


Note. — Information given in this table : (1) ID, (2, 3, 4,5) photo- z, Xred , A v and lcr confidence interval 
from the best SED-fit with a redshift parameter space ranging from z ~0 to 8, (6,7,8) photo-z, Xred and 
A v for the best SED-fit a ssuming a low- 2 : solution (0< 2 : <3), (9, 10) L 1500 deduced from best SED-fit and 
SFR computed using the |Kennicutt et al. | ( fl998| ) relation and corrected for dust extinction. 



Fig. 8. — Distribution of our candidates around the sub-millimeter starburst HFLS3 on a composit image showing the field-of-view 
covered by our HST/WFC3 data. The position of each candidate is displayed by a green circle (red for the candidate in common i n GT C 
and HST samples) and the Herschel starburst by a cyan circle. We also plotted the position of the three faint objects discussed in |4.1| 


However, we noticed that three sources close to HFLS3 
(with 26.5<nipio5VF<28.0) are not detected on ACS data 
(figure [9]). Their break between F814W and F105W is 
<1 mag and their NIR colors could be consistent with 
the Balmer break, therefore we cannot exclude a low -z 
solution for these three objects. Assuming that these 
three sources are at the redshift of HFLS3 , the projected 
distance between them and the Herschel starburst is less 
than 15 kpc. We also compared the overdensity of ob¬ 
jects with brightness similar to these three objects over 


the entire field of view covered by HST. The mean den¬ 
sity over the field is 0.03 object per arcsec 2 whereas it 
is 4x more around HFLS3, reinforcing the overdensity 
hypothesis close to the Herschel starburst. We applied 
the same SED-fitting method described above with the 
same library of templates. We first allowed a large range 
of redshift (0< z <8), but found that all these objects to 
be described SEDs compatible with the SED of a 2^6.3 
galaxy. However, we have to keep in mind that the large 
error bars as well as the small break between F814W 
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and F105W could not exclude a low -z solution for all of 
them. Deep ACS and WFC3 data are needed to strongly 
increase the non-detection constraints in optical and re¬ 
duce the error bars on the WFC3 photometry. 

4.2. Discussion 

In order to constrain the size of a possible overdensity, 
we adopted the galax y fluctuation parameter defined in 


Morselli et al. (2014) 


S=--l 


(9) 


where p is the number of objects selected in our survey 
and p the number of objects expected in a blank field 
covering the same area. 

We then compared the number of objects selected in 
a luminosity range where the completeness is ^100% 
with the number of objects expected in the same inter¬ 
val in magnitude in a blank field. We computed the 
completeness level of the F105W image, where the Zasoo 
is estimated, by adding 10,000 sources per bin of 0.25 
magnitude. We then app lied the extraction parameters 
explained in section |1.3| and compared the number of 
sources extracted with the number of objects added on 
the image. We showed that ^100% of our added ob¬ 
jects are detected up to mpi^w <25.3. Three objects 
in our HST sample are brighter than this magnitude cut. 
We used the UV luminosity functi o n paramete r ization 
published in McLure et al. (|2009 [^Su^t aj^ (2011), 


Bouwens et a T. 


and pouwens et al. ( [2014 ) to es¬ 


timate the number of z~6 galaxies expected m our HST 
survey with mpiosw <25.3 (cf. Table [6]). The expected 
number of objects ranges from 0.5 to TO, therefore the 
over-density parameter ranges from 2.4 to 9.15, with 9.15 
as a strict upper limit on the over-density given the cos¬ 
mic variance and taking into account the contamination 
rate computed in section [2~4| We emphasize that this re¬ 
sult is based on a small field-of-view and a small number 


of objects. Chapman et al. (2009) found a similar value 
(5 ~2.5) for UV-selected galaxies over the GOODS-N 
field, and noticed that this value is lower than that found 
for SMGs in the same field (6 ~10). As shown by Robson 
(2014) there is no over-density of SMGs around 


et al. 


HFLS3. If an over-density occurs around HFLS3 , it is at 
a lower level tha n in the COSMOS AzTEC-3 field (Ca- 
pak et al. |2011|) where they found an over-density factor 
of 11. 

We note that four galaxies among our two samples have 
high SFRs estimated from the SED models (with SFRs 
> 100 M 0 yr -1 ). These sources, however, are und etecte d 
the SC UBA-2 850 pm image of this field (|Ro bson et 
al. 2014). with a noise level of ~ 1.5 mJy/beam. iRiecF 


ers et al. (2014) have studied with ALMA the center of 


the galaxy protocluster associa ted with the z= 5.3 SMG 
AzTEC-3 (Capak et al. ||2011| ), but failed to detect a 
LBG at z= 5.3 with a SFR ~ 20 M 0 yr -1 with ALMA 
data reaching ~ 0.15 mJy/beam (3cr) at 1 mm. Assum¬ 
ing typical Tdust of 35K for the HFLS3 LBGs, a colder 
dust temperature than HFLS3 , observed-frame 850/im 
lies on the Wien part of the SED and the K correction 
from z=5.3 to z=6 .3 remains positive. Scaling the SED 
and depths reached between ALMA at 2 : = 5.3 at 1 mm 
and SCUBA-2 for galaxies at 2 ^ 6.3 at 850 /rm, we 


find that the SCUBA-2 depth is only adequate to detect 
galaxies above 3 a with instantaneous SFRs > 500 M 0 
yr -1 . This ignores the cosmic microwave background 
(CMB) that will absorb a few per cent of the sub-mm 
flux at z=6 .3 compared to 2:=5.3. Therefore SCUBA-2 


observations of the depth of the Robson et al. (2014) 


data are unlikely to detect galaxies with similar SFRs 
and SED properties as LBGs in the AzTEC-3 protoclus¬ 
ter. 

4.3. Distribution in luminosity of the z~6 candidates 

The redshift interval covered by our survey is sufficient 
to compute the UV luminosity function at z~6. This 
again allows a way to study if there is an overdensity of 
objects at a given luminosity. Several methods have been 
devel oped to compute the UV LF at different red shifts 
Bouwens et al. 2012b, Willott et al. 2013). We 


made use of the previous SED-fitting work to compute 
the number density of objects at z~6 as a function of the 
UV luminosity (see |Laporte et al.~^|2015 for details). 

The effective surface explored by our HST survey is 
computed by masking the bright objects on the detection 
image. The comoving volume explored between z ^5.5 
and 2^6.5 is V ~10 980 Mpc 3 . The number densities 
we found are shown on Table [5] and the error bars take 
into account the cosmic variance. These numbers are 
consistent wit h previous findings in this redshift interval 
(e.g. Willott et al. 2013 


Bouwens et al. 


2003. 


We adopte d the Schechter parameterization (|$chechter 


et al. 1976) of the UV LF defined by: 


<f>(M) = 


MM (1 „ -0.4(M-M*))“+1 exp 01 q-O.4(M-M*)) 

2.5 

( 10 ) 

where M*, <F* and a are the 3 Schechter parameters to 
be adjusted. 

In order to show the influence of the densities com¬ 
puted in this study, we fitted the shape of th e UV LF 
by co mbining them with the density from Bouwens et 


al. (2014) in the faint-end slope and in the bright 


end with those published by Willott et al. (120131). 
We used a y 2 minimization with parameters ranging 
from 10 -5 to 10 -2 Mpc -3 .mag -1 , -22.0 to -19.0 and - 
2.2 to -1.1 respectively for <F*, M* and a. The follow¬ 
ing parameterization is found at z~6 : M*=-20.17 


■+ 0.35 

- 0.155 


4V=(L21I 019 )xlO -3 Mpc -3 .mag -1 and <a=-1.77±0.23. 
Error bars on each parameters are deduced from the la 
confidence interval. These parame ters are in good agree¬ 
ment with prev ious results (e.g. [B ouwens et al. ||2012| 
Su et al. 1|2011 and McLure et al. |[20Q? - Table [fif. The 
shape of the UV LF deduced from this study as well as 
the number densities we computed are plotted on Figure 
PH We find no evidence for an over-density of galaxies, 
above the field LF, at any of the luminosities probed by 
the data. 

5. CONCLUSIONS 

We have presented in this paper the results of the pho¬ 
tometric analysis of the environment of a sub-millimeter 
starburst at z~6. 34, combining both wide area ground- 
based data and high-level quality HST data to explore 
a large range of luminosities. We applied the Lyman- 
break galaxy selection technique and found 10 galax¬ 
ies that are at z~6. This sample includes two from 
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Fig. 9. — Thu mbnail images of the faint sources located around HFLS3. The size of each stamp is ~8.3” x8.3”. The lens system studied 
in |Cooray et a l. ([2014) is at the center of the field. The position of each faint candidate is display by a circle. 
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Fig. 10. — Result s of the Voronoi tessellation analysis following the method described in Ramella et al. | (2001) on the sample of faint 
sources described in |4.1| The background threshold above which fluctuations are consistent wit! an over-dense held is estimated at 1.74 
per Voronoi cell, which is not a significant detection of an over-density in the imaging data. 


our ground-based search (m z ~ 25) and eight more from 
our search with HST/WFC3 data (with mpiosw h 1 the 
range 24.5 to 26.5). The size of each sample is well- 
consistent with expectations from previous findings using 
blank field surveys and seems incompatible with an over¬ 
density of luminous (rriFios <25.9) galaxies. We used a 
standard method to estimate the photometric proper¬ 
ties of each source and we used these SED-fitting results 
to compute the UV LF at z~6. The parameterization 


of the Schechter function we deduced (M*=-20.17^0.15, 

4>*=(1.2llong) x 10 -3 Mpc -3 and <a=-1.77±0.23) is in 
good agreement with previous finding in this redshift in¬ 
terval. 

We do not find any strong evidence for HFLS3 being a 
member of a proto-cluster of luminous galaxies as is the 
case of the well known SMG Aztec-3 (Capak et al. 2011). 
There the overdensity parameter was found to be about 
11 ; in the case of HFLS3 , we place an upper limit on the 





















12 


Laporte et al. 



Fig. 1 1.— Shape of the UV luminosity function a t z ~6. Number densities from this study are plotted in red, previous finding are also 
plotted (|McLure et al. ||2009| |Bouwens et al. ||2007| and |Willott et al. |2013|. The parameterization computed from the study of HFLS3 
envir onment is displayed by the black line, the ot her l ines display the parameterization published in previous references (dashed line from 
|Bouwens et al. |2012||and dotted-dashed line from|McLure et al. [2009). The la confidence intervals on each parameter are also shown on 
the lower panels. 


TABLE 5 

Number density at z~6 


A/1500 

<F(Mi50o) 

d<F 


xl0 -4 [Mpc -3 .mag -1 ] 

X 10 _4 [Mpc _3 .mag _1 ] 

-21.75T0.250 

0.09 

+0.23 

-0.09 

-21.25±0.250 

1.01 

i0.49 

-20.75i0.500 

2.60 

il.22 


Note. — Number density computed following a method using the 
redshift probability distribution. Error bars include d poisson uncer- 
tainties and cosmic variance computed from |Trenti 8z Stiavelli ( 2008]) 


overdensity of ~9 after taking into account the cosmic 
variance of existing z ~ 6 LF measurements. The lack 
of a significat overdensity is also confirmed by a Voronoi 
tessellation analysis that included all the faint objects 
fulfilling the color-criteria defined for z~6 objects, but 
without a well-defined break between the optical ACS 
and near-IR WFC3 imaging data. 

However, we noticed at least three faint objects within 
3” from HFLS3. They are undetected in ACS images 
and are detected on WFC3 data. If the redshift of these 
sources are confirmed at then the Herschel star- 

burst is located in an over-dense region composed by faint 
objects but with an extent of 36 kpc. Deeper data com¬ 
bined to spectroscopic observations are needed to asses 
these conclusions. It is more likely that these are associ¬ 
ations with the merger system that might be triggering 
the starburst. 
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TABLE 6 

UV LF PARAMETERIZATION AT 



Reference 

M* 

4>* 

X 10 -3 [Mpc _3 .mag _1 

a. 

] 


This work 

-20-16t°;?j| 

1 36+ 0 - 40 
l.OO — o i9 

-1.70T0.23 

iMcLure et al. 

(2009) 

-20.04T0.12 

1.80T0.50 

-l.71T0.il 


Su et al. 

(2011) 

-20.25T0.23 

i yy+0.62 

i ‘ ' 1 -0.49 

-1.87T0.14 

Bo 

uwens et a 

1. 

(2012 

-20.37T0.30 

1 4 +1-1 
-*-•^-0.6 

-1.73zb0.22 

Bouwens et al. 

(2014 

-20.93T0.25 

0 49+ 0 - 26 
u -^ y -0.17 

-1.85T0.10 








Note. — Parameterization of the Schechter function presented in 

this paper and 


published by other teams. 





Fig. 12 .— SED-fitting of the 3 faint objects highlighted around HFLS3 using two assumptions on the redshift: (black line) allowing 
a large range of redshift 0< z <8 ( X^ e d 0.02, 0.01 respectively for ID1, ID2 and ID3) and (orange line) fixing the redshift at the 

redshift of HFLS3 (Xred M).15,0.06, 0.01 respectively for ID1, ID2 and ID3). The blue dots show the photometry of these 3 faint objects 
(SExtractor MAG_AUTO) and the 2cr upper limits in case of non-detection. Error bars are computed using the noise measured in empty 
apertures around the objects. 


age.This work is based on observations made with the 
Gran Telescopio Canarias (GTC), instaled in the Span¬ 
ish Observatorio del Roque de los Muchachos of the In¬ 
stitute de Astroffsica de Canarias, in the island of La 
Palma. The William Herschel Telescope is operated on 
the island of La Palma by the Isaac Newton Group in the 


Spanish Observatorio del Roque de los Muchachos of the 
Institute de Astrofsica de Canarias.The data presented in 
this paper will be released through the Herschel Database 
in Marseille HeDaM (hedam.oamp.fr/HerMES) 
Facilities: GTC (OSIRIS), HST (ACS), HST (WFC3), 
Spitzer (IRAC), WHT (LIRIS). 
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